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!~r: Abstract 

P^ ' It is believed that choked gamma-ray bursts (CGRBs) are the potential candidates for the 



Oh! 
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production of high energy neutrinos in GeV-TeV energy range. These CGRBs out number the 



O . successful GRBs by many orders. So it is important to observe neutrinos from these cosmological 

^ ' objects with the presently operating neutrino telescope IceCube. We study the three flavor neutrino 

oscillation of these high energy neutrinos in the presupernova star environment which is responsible 

K* . for the CGRB. For the presupernova star we consider three different models and calculate the 

^^ ' neutrino oscillation probabilities, as well as neutrino flux on the surface of these star. The matter 

'^ ' effect modifies the neutrino flux of different flavors on the surface of the star. We have also 

^D ' calculated the flux of these high energy neutrinos on the surface of the Earth. We found that 

m 

for neutrino energies below < 10 TeV the flux ratio does not amount to 1:1:1, whereas for higher 
energy neutrinos it does. 



PACS numbers: 14.60.Pq; 98.70.Rz 



massive stars are related to long type 
5|, |6| . The core collapse of massive stars 



I. INTRODUCTION 

The long Type (> 2 s) of gamma-ray bursts (GRBs) constitute about 3/4 of the total 
observed GRBs. The observed correlations of the following GRBs with supernovae (SNe), 
GRB980425/SN 1998bw [1|,GRB 030329/SN 2003dh [2], GRB 031203/SN 20031w[3] and 
GRB 021211/SN 20021t[4 show that core collapse o 
of GRBs. These mostly occur in star forming regions 
resulting in a relativistic jet which breaks through the stellar envelope is a widely discussed 
scenario for gamma-ray burst production. In this scenario, the gamma-rays are produced by 
synchrotron and/or inverse-Compton scattering of Fermi accelerated electrons in optically 
thin shock, when the jet has emerged out of the stellar envelope. These same shocks are 
also responsible for the acceleration of the protons into relativistic velocities and collisions 

Q 

of these with the MeV photons produce neutrinos of energy ~ 100 TeV[7|]. 

The formation of jets may be a common phenomena in a collapsar scenario and depending 
on the composition (baryon load) all the jets may not have sufficient energy and momentum 
to punch trough the stellar envelope. From the observed rate of GRBs, only about < 10~^ 
core collapse supernovae produce highly relativistic jets (Lorentz factor F > 100) which can 
penetrate through the envelope and produce GRB events in an optically thin environment 
outside the star 8}]. On the other hand, a large fraction of them will fail to emerge out 
of the envelope, which may give rise to orphan radio afterglow instead of 7-ray emission. 
However when the jets are making their way through the star they can accelerate protons to 
energy > 10^ GeV through internal shocks well inside the stellar envelope. Also the buried 
jet produces thermal X-ray at ~ 1 keV which acts as the target for the delta-resonance 
to produce > 5 TeV energy neutrinos through photopion production, which penetrates 
through the envelope(9(]. There can also be neutrino production due to pp and pn collisions 
involving relativistic protons from the buried jet and the thermal nucleons from the jet 



and the surrounding 10| . However, when the jets are still making their way through the 
star, precursor TeV neutrinos are inherent to the collapse of a massive star irrespective of 
whether GRB is produced or not. Forward moving jets which are unable to emerge out of 
the envelope may be rich in baryons which make them mildly relativistic, but at the same 
time efficient producers of TeV neutrinos through photomeson interaction. The detection 



of low luminosity GRB 060218 



ll| suggests that their number (low luminosity GRBs and 



dark GRBs) may be quite large compared to the hi gh l uminosity GRBs, which can probably 
contribute more to the TeV neutrino background 12| than the high luminosity ones and 



could be detectable by present day neutrino telescopes (e.g. in IceCube) which can shed 
more light on the nature of the central engine, as well as the acceleration mechanism of 
high e,...g,v oos,nic .ays i. the p.esupen>ova stafiQ. FVo,„ a. ,„d,v,dua. co„ap.e/GRB 
burst at a distance z ~ 1, about 0.1-10 upward going muon events can be detected in a km^ 
detector [qI]. 

These high energy neutrinos propagating through the presupernova star with a heavy 
envelop can oscillate to other flavors due to the matter effect. In fact the mater effect is 
well known for neutrinos propagating in the sun, supernova as well as in the early Universe. 
In optically thick hidden sources where gamma-rays are not observed directly and ~ TeV 
neutrinos are produced due to p7, pp and pn collision, the flux ratio at the production site 
and on the surface of the star may be different due to matter effect on their oscillation. 
Recently for two neutrino flavors it is shown that for the choked GRBs, the multi-TeV 
neutrino signals proposed by Meszaros and Waxman |9J can undergo substantial resonance 
oscillation before escaping from the He envelope if the neutrino oscillation parameters are 
in the atmospheric neutrino oscillation range. This would alter the neutrino flavor ratio 
escaping from the stellar envelope, and subsequently the detected flavor ratio on Earth 15 1. 
So in this context it is important to study the matter effect of the presupernova star on the 
oscillation of high energy neutrinos emerging out of it. 

The neutrino oscillation in vacuum and matter has been discussed extensively for solar, 
atmospheric as well as accelerator and reactor experiments. Models of three flavor neutrino 
oscillations in constant matter density 16Nl8| , linearl y v arying density 19|, |20| and exponen- 
tially varying density 211] have been studied. In Ref. 22| T. Ohlsson and H. Snellman have 
developed an analytic formalism for the oscillation of three flavor neutrinos in the matter 
background with varying density, where they use the plane wave approximation for the neu- 
trinos (henceforth we refer to this as OS formalism). Here the evolution operator and the 
transition probabilities are expressed as functions of the vacuum mass square differences, 
vacuum mixing angles and the matter density parameter. As application of the above for- 
malism, the authors have studied the neutrino oscillations traver sing the Earth and the Sun 
for constant, step-function and varying matter density profiles 23|,|2J|. To handel the varying 
density, the distance is divided into equidistance slices and in each slice the matter density 



is assumed to be constant. In these calculations they have considered the CP phase to be 
real by taking the phase factor 5cp = so that the neutrino mixing matrix is real. 

Although the OS formalism is simple and used for low energy neutrino oscillation, so far 
it has not been used to study the propagation of high energy neutrinos neither in the stellar 
envelope where the density is high nor in the Earth. In the present work we are using this 
formalism to study the three flavor high energy (energy in the range 100 GeV to 100 TeV 
) neutrino oscillation when traversing the presupernova star medium and reaching to the 
Earth by undergoing vacuum oscillation in the intergalactic medium. 

The paper is organized as follow: In SecUl] we review the OS formalism used for the 
calculation of neutrino oscillation probabilities in a medium. In Sec. Illll we discuss about the 
presupernova star models to explain the choked GRBs. The discussion of our results are 
given in Sec. IIVI followed by a summary in Sec. |Vl 

II. NEUTRINO OSCILLATION FORMALISM 

In this section we shall summarize the formalism used by OS [22] for the calculation of 
the oscillation probability of the three active flavors. A flavor neutrino state is a linear 
superposition of mass eigenstate and is given as 



) = Ef^-i^«)' (1) 



a=l 



where a = e^n^r (flavor eigenstates) and a = 1,2,3 (mass eigenstates). The matrix Uaa, is 
the three neutrino mixing matrix given by. 
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Uel Ue2 Ue3 

Uf^l U^2 f^^3 

\Url Ur2 Ur3l 



' C13C12 C13S12 sise '^'^p^ 

"■S12C23 ~ Ci2S23Si3e* "P C12C23 — Si2S23Si3e* "'^ >S23Ci3 
S23S12 - C23Si3Ci2e'^'"' -S23C12 - Si3Si2C23e*'''='' C23C13 



(2) 



y 623S12 - C23Sl3Cl2e-'"' -S23C12 - 613612^236 '-" C23C13 y 

where Cij = cos 6ij and Sij = sin 6ij for i,j = 1,2,3. The neutrino mixing angles are ^12,^13 
and ^23- The 6cp is the CP violating phase. As the CP violation is not observed in the 
neutrino sector, we put this phase 6cp = in our calculation. 

While the neutrinos travel from the production point to the detection point, the flavor 
ratios will evolve as a result of their oscillations. These neutrinos will go through vac- 
uum and matter during their propagation. In vacuum, the Hamiltonian that described the 



v 



propagation of the neutrinos in the mass eigenstate basis is described by 

^ El Q 0^ 
Hm= E2 , (3) 

Ey 

where Ei, for i = 1,2,3 refer to the energy of each neutrino mass eigenstate. This Hamilto- 
nian can be written in the flavor basis through the unitary transformation described by the 
matrix U from equation ([2]), as 

Hf = UH^U-\ (4) 

In matter, the Ug and Ue will interact through both charge and neutral current, whereas 
z/^ and Ur and their anti-neutrinos will interact through neutral current only. So for the 
oscillation of electron neutrino (or anti-neutrino) to other flavor, only the charge current 
term will contribute. For the oscillation of z/^ ■H- Vri there is no contribution from the 
matter up to leading order. Thus the matter contribution to the Hamiltonian in the flavor 
basis can be expressed as 



^A 0^ 



Vf 








(5) 



where A represents the potential due to the interaction of Ve/^e with matter and is given by 



A = ±V2GFNe = ±V2G 



P 



(6) 



where Tn^v is the nucleon mass, Gp is the Fermi coupling constant and p is the matter density. 
The neutral current contribution to the potential is the same for all the three neutrinos, so 
here we do not take that into account. In the mass basis, the total Hamiltonian is given by 

Hm = Hra + U~'VfU (7) 

The total Hamiltonian in the flavor basis is written as 



Uf = Hj + Vj. 



(8) 



For neutrino propagation in a medium, the Hamiltonian is not diagonal, neither in the mass 
basis nor in the flavor basis, so one has to calculate the evolution operator in any of these 
basis. 



In the mass basis, the evolution of the state at a later time t will be obtained by solving 
the Schrodringer equation 

^^^=^„.ka(t)), (9) 

and the solution to this equation can be expressed in terms of the evolution operator as 

ka(t))=e-«-*|z/a(0)) (10) 

= f/„(t)|z/,(0)), 

where Um{t) is the evolution operator in the mass basis and in the flavor basis this can be 
written as 

Uf{t) = uuut)u-\ (11) 

As neutrinos are relativistic, we can replace t by the path length L, where we use the natural 
units c = 1 and h = 1 . To obtain Um{L) we have to evaluate the exponential of the matrix 
—iTimL and also introduce a traceless matrix T deflned as 

T = «,..-^!l^^ (12) 

The trace of the Hamiltonian in the mass basis is 

3 

tr(?/J = ^E, + A (13) 

i 

For the square matrix, T, its exponential is given by 

°^ rpn 



n=0 



Using Cayley-Hamilton's theorem, this inflnite sum can be expressed as a flnite sum and 
can be given by 



e 



n=0 

where A^ is the dimension of the matrix T. Here we have A^ = 3 and the matrix T has three 
eigenvalues Aj with i = 1,2,3. The solution of the characteristic equation of the T matrix 
will be obtained by solving 

A=^ + csA^ + ciA + Co = 0. (16) 

The coefficients of A are given as 

C2 = -tr{T) =0, Co = -det(T), (17) 



and 

Ci = T11T22 — T12 + T11T33 — T^3 + T22T33 — T23. (18) 

The eigenvalues are given by 



Ai 



21/332/3 X 



2.3 22/331/3X 2 X 21/332/3 ' V ) 



(20) 



with 

X= fy3y^4cf + 27c2-9coj . 

With the use of the above, the evolution operator in the mass basis can be written as 

f/„(L) = e'^^-^ (21) 

_ . ^^ -.Lx. [(Ag + ci)/ + A„T + T2] 
^^, 3A2 + C1 

a=l " 

where = e-^L*'^i'Hm)i/3 jg g^ complex phase factor and / is the identity matrix. The evolution 
operator in the flavor basis is given by 



Uf{L) = e-*^^^ (22) 

(A2 + ci)/+A,T + T2l 



3 

-iLXa 



■Y. 



e 



3\l + ci 

a=l " 

where T = UTU^^ and T^ = UT'^U^^ . The probability of flavor change from a to /3 due to 
neutrino oscillation through a distance L can be given by 



P,^^,^{L)^P^p{L) = \{vp\Uf{L)\v^)\ 



2 



3 3 



^aP - 4^ '^Pa{L)paPb{L)pa siu^ Xab, (23) 



a=l 6=1 
a<b 



where we have defined 



D/r\ (A^ + Ci)5/3a + AaT^a +T|^ 
i^a(^)/3« - ^^^^^ . (24) 



The matrices Tp^ and T'o^ are symmetric and defined as 

3 3 
Tap = Tj3a = 2_^ Z-^ UaaUj3bTab, (25) 



a=l fe=l 



and 



rp2 

-'-a/3 



rp2 

^ 13a 



Yl E UaaU^hT, 



•2 
ab' 



Also we have defined the quantity 



•^ab 



a=l 6=1 



(Aa - Xb)L 



(26) 



(27) 



The matrix T is written explicitly as 

/ 



T, 



ab 



Tu AUelUe2 AUelUeS 

AUelUe2 T22 AUe2Ue3 

AUelU,, AU,2Ue3 T33 



(28) 



where the diagonal elements of the above matrix are given by 

3 



Taa = AUl + hY^ Eab-A). 



b^a=l 



Here E, 



ab 



-E, 



ba 



Ea — Eh and the energies satisfy the relation 
E12 + E23 + E31 = 0. 
The neutrino oscillation probabilities satisfy the condition 

y^ Pai3 = 1, for a,(3 = e, fi, r, 

and a similar condition is satisfied for anti-neutrinos which we define as P; 



(29) 



(30) 



(31) 



Q/3- 



Using the Eqs. (!23l) and ( l3Tl) we can calculate the probability of transition from one flavor 
to another. By substituting A = in Eq.(l23l) we get the vacuum probability. For matter 
with varying density the distance L can be discretized into small segments with a constant 
density in each segment and use this formulation repeatedly in each segment. By doing so we 
can study numerically the neutrino oscillation in any type of density profile. For neutrinos 
traversing a series of matter densities pi ioi i = 1 to n, with their corresponding thickness 
Li, the total evolution operator is the ordered product and is given as 



Uf{L) = llUf{L, 



(32) 



where Y17 Li = L. For rapidly changing profiles one has to consider smaller segments to 
make sure that density is almost constant in each of the segments. In a series of papers by 
OS, this method has been applied for different density profiles of the Sun and the Earth, to 
study the MeV energy neutrino oscillation. 
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III. PRESUPERNOVA STAR MODELS 



The popular models for the long duration gamma-ray bursts (LGRBs) are the core col- 



lapse supernova models called the collapsar models, where the core of a massive star co 



lapses 



to form a black hole and drive an ultra relativistic jet, which breaks out of the star[25|]. Due 
to the rotation of the star, the mass distribution along the rotation axis is lower than the 
equatorial region which helps to launch the relativistic jet along the rotation axis of the pro- 
genitor forming a black hole at the center. The jet is probably powered by the annihilation 
of neutrino anti-neutrino pairs or some other electromagnetic process. Even though the for- 
mation of jets may be a common phenomena in a collapsar scenario, the mass density of the 
progenitor is quite high, so it is not so trivial that all the formed outflow can always be col- 
limated and punch through the stellar envelope as a jet. Depending on the initial mass and 
metallicity, the presupernova star can have different compositions and different radii which 
can not be probed observationally. So for this reason we resort to models that assume some 
standard processes considered by different authors, and which seem to fit the observational 
data well. A general model for the density of the progenitor of the SN is based on a power 
law, p oc r~"[26] and for a star this parametrization is valid only in a range of r, but should 
be a fairly good approximation for our calculation. The power index n depends on the stage 
in wich the progenitor is. There are two commonly used values: n = 3 and 17/7. The value 



n = 3 
model 



27 



resulting from the study of the SN1987A and 17/7 for a simple blast wave 
29|. Also depending on the composition of the presupernova star, the spectra will 
have some specific features that allow a classification of the SNe. The presupernova could 
be a red super giant (RSG) or a blue super giant (BSG), both having hydrogen envelope 
but different energy transfer mechanisms (convective or radiative respectively). 

Based on above and with the fact that some of the Long GRBs observed are related 



to SNe (Type Ib/c), we are going to use three possible models as considered in [30|. The 
characteristics of these progenitor models are that they are having an iron core of radius 
rpe ~ 10^ cm surrounded by a He core extending up to thc ~ 10^^ cm where the density 
is pHe ~ 10~^ g/cm^. In some cases a hydrogen envelope surrounds the He core extending 
to vh ^ 10^^ cm with a density of pn ~ lO""^ g/cm^ . The presupernova stars which are 
believed to be the strong contender for the Long GRBs are Type Ic SNe which have lost 
their hydrogen envelopes as well as most of the He envelope before the explosion. These 




r (cm) 



FIG. 1: Density profiles of tlie pressupernova star ( a blue supergiant) of models A, B and 
C with a radius R^, = 3 x 10^^ cm. In these models, the high energy neutrinos are 

produced at a radius Vj = 10^°'^ cm. 

objects are not interesting from the point of view of neutrino oscillation because their radii 
are too small to have any appreciable effect. On the other hand, the presupernova stars 
with the He envelope (Type lb SNe) and even the H envelope (Type II SNe) intact will be 
favorable for TeV neutrino production as well as their oscillation in the stellar environment. 
We consider three different presupernova models as shown in FIG.[T1 In all these models the 



radius of the BSG is taken to be i?^, ~ 3 x 10 cm. The jet evolves at a radius r 



10 



10. g 



cm < rHe and also the > 5 TeV neutrinos are produced at a point between rj and rne so 
that the jet can acquire relativistic velocity on the surface of the He envelope. The models 
are: 



Model A 



pir) = p, X ( ^ - 1 



cm-" 



(33) 



This model represents a star with a radiative envelope. It has a polytropic structure 
with a polytropic index n = 3 and the characteristic density p* = 4.0 x 10~^ g/cm^. 
This model is valid only in the region of the star that lies between the point were the 
jet is produced and the star envelope i.e. r^- < r < i?*. 
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Model B 

' [^Y ;10i°-Scm<r <n g 



p{r) = p,x 

I fil^V 

(rt-fl.)' 

(34) 



I r J (r.-R.)" ' ^ > ^fe 



This model is for a BSG with polytropic index n=17/7 and with r^ = 10^^ cm. The 
density has a power law behavior with a characteristic density p^, = 3.4 x 10^^ g/cm^. 
This profile was obtained from the fit to SN1987A data. 



Model C 



p{r)=p,A{^-l] ^, (35) 



r / cm^ 



where the parameters of model C are given as 

(2.1, 20) ;10^°-*cm<r < 10^^ 

{neff,A) = < 

(2.5, 1.0) ;r > 10^^ cm. 



(36) 



This model has two free parameters ng// and A that can be fitted to produce a drop in 
the density after the hehum core. While the parameter A can be used to set a density 



drop at the edge of the helium core, rig// gives the effective polytropic index 26|]. The 
characteristic density here is p* = 6.3 x 10~^ g/cm^. The density profiles of all these 
three models are shown in FIG. [1] 

IV. RESULTS 

We use the neutrino oscillation formalism of OS given in Sec. II and calculate the neu- 
trino oscillation probabilities and fluxes in presupernova star models A, B and C and also 
on the surface of the Earth after the neutrinos have undergone vacuum oscillation in the 
intergalactic medium. For all these calculations we take the neutrino energy in the range 
100 GeV to 100 TeV , although it may be difficult to produce neutrinos with energy more 
than 10-20 TeV in the presupernova star environment. 

11 







FIG. 2: The oscillation probabilities for neutrinos (black solid line) and anti-neutrinos (red 

dotted line) are plotted as functions of neutrino/anti-neutrino energy E^, for model A, 

Set-I where first figure is for Pee and Pgg, second figure is for P^^ and Pep,, third figure is for 

Per, Pef and the fourth one is for P^,- and Pp,f respectively. 



We use the standard neutrino oscillation parameters obtained from different experiments 
for analysis of our results. The neutrino parameters used are as follows: 



Ami = 8.0 X lO^^el/l 6,0 = 33.8° and 6 



''21 



''23 



45° 



Am 



32 



3.2 X 10"^ ey^ ^13 = 8.8° and 6cp = 0. 



(37) 



We also we take ^13 = 12° to observe the variation in oscillation probability due to change 
in this angle. 

For the numerical calculation, we divide the distance (P* — r^) into small slices, each with 
a constant density and calculate Uf{L) for each individual slices as discussed in Eq.f p2|) . 
Afterward we use Eq. (!23l) to calculate the probabilities. We let the high energy neutrinos 
propagate from the production point at Tj = 10^°'^ cm towards the surface of the star 
P*, where we calculate their survival and transition probabilities Pa/3 and P^^, as well as 
fluxes. We consider two different sets of parameters: P* = 3 x 10^^ cm and 613 = 8° (Set-I); 
P^, = 2.7 X 10^^ cm and ^13 = 12° (Set-II); to observe the variation in the probabilities at 
different depth from the star surface and different ^13. 

The survival and transition probabilities of u and z/ for the parameter Set-I, for models 
A, B and C are shown in FIGs. [21 [3] and H] respectively. We have also plotted the variation 

12 




FIG. 3: This is same as FIGj2]but for model-B. (See the electronic edition of the Journal 
for a color version of the figures and the color specification is same as FIG 12]). 




FIG. 4: This is same as FIG J2] but for model-C. 



in oscillation probabilities for Set-II in FIG. [51 We observe that for a given neutrino energy 
El,, Pgg is always above Pee and the probabilities are highly oscillatory for E^ < 10^^ eV. 
For El, ^2 X 10^^ eV, both Pee and Peg increase towards unity. While the increase in Pgg is 
smooth, the increase in Pee is accompanied by a rapid oscillation, as shown in these figures. 
By reducing the P* to 2.7 x 10^^ cm and 6*13 = 12° i.e. Set-II (FIG. El), we observe that, 
although there is variation of probability on the surface, the overall behavior is exactly the 
same as the parameter Set-I shown in FIG 121 So for our further discussion we only consider 
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FIG. 5: This is same as FIGj2]but for the parameter set-II. 
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FIG. 6: The energy difference is plotted as a function of potential energy A to look for the 

existence of resonance. 

the parameter Set-I for the analysis of our results. All the transition probabilities P^.^^ Per, 
PfiT Pep., Per and P/if are highly oscillatory in all the energy ranges (mostly < 10^^ eV). 
For neutrino energy E,^ > 3 x 10^^ eV the transition probabilities Pe^, Per, Pep. and Per 
go to zero, which shows that the medium has almost no effect on high energy neutrinos. 
On the other hand P^r and P^^ are highly oscillatory in this energy range. The Pg/^ and 
Per are different from each other due to the matter effect. Comparison of Pgp and Pgf 
shows that they oscillate with same amplitude, but are out of phase by 180°. The transition 
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FIG. 7: Neutrino and anti-neutrino fluxes on the surface of the star. In (a), (b) and (c) 

black solid line is for ^^^, blue dashed lines is for ^^^, and red dotted line is for $,^^. In (d) 

black solid curve, blue dashed curve and red dotted curve are for $p^, $p^ and $p^ 

respectively. 



(a) Vacuum (inside star) 
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FIG. 8: Same as FIG|7]but the flux is calculated on the surface of the Earth. 



probabilities P^r and Ppf are different for energies below ~ 3 x 10 eV, but are almost the 
same above this energy for all the models A, B, and C (last of FIGs. El |3], [Hand [5]) which 
is due to the negligible medium effect on the high energy neutrinos. 

The energy eigenvalues of the neutrinos in the matter are A^ and the energy difference 
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FIG. 9: The total flux (neutrino+anti-neutrino) on the surface of the star, where the black 
solid, blue dashed and red dotted curves are for $j,^, $i^ and $j,^ respectively. 
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FIG. 10: The total flux (neutrino+anti-neutrino) on the surface of the Earth, where the 
black solid, blue dashed and red dotted curves are for $i,^, ^^ and ^^^ respectively. 



|Aa — Afol is related to the effective mass square difference 



An — A;, 






(38) 



In FIGJSl fo^^ the illustrative purpose, we have taken E^ = 1 TeV with parameter Set-I 
and model-A to show the resonance position as a function of matter potential A. It shows 
that there is only one resonance position around yl ^ 2 x 10~^^ eV (between |Ai — A2I and 
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1^1 — ^sD- By taking E^, = 10 TeV we found that the resonance position is almost the same. 
The anti-neutrinos will not satisfy the resonance condition because of the change in the sign 
of the potential A in Xa- 

In the mildly relativistic jet internal shocks can develop and accelerate protons to very 
high energies. These protons would interact with the ~ keV thermal X-ray photons to 
produce TeV neutrinos via the process p + 7 — !■ A+ -^n + TT^^n + fi^ + i'fj,-^n + e~^ + 
v^i + Ve + ^11- In the above process the standard neutrino flux ratio at the production point 
is $^^ : $° : $°^ = 1:2:0 ($^^ corresponds to the sum of neutrino and anti-neutrino flux 
at the source). The flux observed at a distance is given by 

^u^ = Y. ^%P<-P^ a, /3 = e, /i, r. (39) 

Using models A, B and C, we have calculated the normalized fluxes of neutrinos and 
anti-neutrinos on the surface of the presupernova star as well as on the surface of the Earth, 
which are shown in FIGs. [7] and [HI In these flgures, (a) is the flux calculation in vacuum 
(where the matter potential A = Q inside the star) and the (b), (c) are for models A, B 
and (d) is for anti-neutrino fluxes for model A respectively. The comparison of the vacuum 
oscillation i.e. (a) with the rest i.e. (b), (c) and (d), shows that for E'j, > 6 x 10^^ eV both 
the fluxes are the same, which signifles that above this energy the matter effect is negligible 
and below this energy the normalized flux is oscillatory and energy dependent. Also the 
comparison of fluxes of neutrinos (b) and anti-neutrinos (d) in model A below -Ej^ > 6 x 10^^ 
eV are different. 

In FIGs. [9] and [To] we have shown the total flux (neutrino and anti- neutrino) of z/g, z/^ and 
Vj. on the surface of the star and on Earth. It can be observed that the fluxes of different 
neutrinos are different on the surface of the star. For E'^ < 4 x 10^^ eV, the fluxes on Earth 
are also different but above this energy the flux ratio is close to 1:1:1. 

V. SUMMARY 

It is observed that only a very small fraction of core collapse SNe are responsible for 
GRBs and majority of them fail to produce 7-rays. This majority class are efficient emitters 
of TeV neutrinos which can be observable in IceCube. Here we study the matter effect of the 
presupernova star on the propagation of these high energy neutrinos by taking into account 
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all the three active flavors and using the formalism developed by OS. We observed that the 
neutrino oscillation within the presupernova star depends on the neutrino energy and alter 
their flux on the surface of the star. We also calculated the fluxes of these neutrinos on 
the Earth after they travelled the intergalactic medium (vacuum oscillation) and found that 
low energy neutrinos E,^ < A TeV have different fluxes whereas above this energy the flux 
ration is close to 1:1:1. So, possible detection of high energy neutrinos with E^ < 10 TeV by 
IceCube can be important due to the change in their fluxes when they propagate through 
the choked environment of the presupernova stars. Also detection of these neutrinos might 
shed more light on the type of progenitors and the acceleration mechanisms of the high 
energy cosmic rays. As application and continuation of the present work, the calculation of 
the neutrino induced muon (track) to electron (shower) ratio in IceCube is in preparation. 
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bia, for their kind hospitality during his several visits. We thank Karla Patricia Varela for 
helpful discussions. This work is partially supported by DGAPA-UNAM (Mexico) Project 
No. IN103812. 



[1] T. J. Galaina, P. M. Vreeswijk, J. van Paradijs, C. Kouveliotou, T. Augusteijn, O. R. Hainaut, 



F. Patat and H. Bohnhardt et al, Nature 395, 670 (1998) (astro-ph/9806175l . 
[2] P. A. Mazzali, J. S. Deng, N. Tominaga, K. Maeda, K. .Nonioto, T. Matheson, K. S. Kawabata 

and K. Z. Stanek et al, Astrophys. J. 599, L95 (2003) |astro-ph/0309555]. 
[3] B. Thomsen, J. Hjorth, D. Watson, J. Gorosabel, J. P. U. Fynbo, B. L. Jensen, M. I. Andersen 



and T. H. Ball et al, Astron. Astrophys. 419, L21 (2004) |a stro-ph/040345l 1. 
[4] M. Delia Valle, D. Malesani, S. Benetti, V. Testa, M. Hamuy, L. A. Antonelli, G. Chincarini 



and G. Cocozza et al, Astron. Astrophys. 406, L33 (2003) | astro-ph/0306298| . 
[5] S. E. Woosley and J. S. Bloom, Ann. Rev. Astron. Astrophys. 44, 507 (2006) 

|astro-ph/060914 2]. 
[6] B. Zhang, B. B. Zhang, F. J. Virgih, E. W. Liang, D. A. Kann, X. F. Wu, D. Proga and 



H. J. Lv et al, Astrophys. J. 703, 1696 (2009) [arXiv:0902.2419;' [astro-ph.HE]] . 
[7] E. Waxman and J. N. Bahcall, Phys. Rev. Lett. 78, 2292 (1997); E. Waxman and J. N. Bahcall, 
Phys. Rev. D 59, 023002 (1999). 

18 



[8] E. Berger, S. R. Kulkarni, D. A. Frail and A. M. Soderberg, Astrophys. J. 599, 408 (2003); 
[9] P. Meszaros and E. Waxman, Phys. Rev. Lett. 87, 171102 (2001). 
[10] S. Razzaque, P. Meszaros and E. Waxman, Phys. Rev. Lett. 93, 181101 (2004) [Erratum-ibid. 

94, 109903 (2005)]. 
[11] T. J. Galama et al, Nature, 395, 670 (1998); S. Campana et al, Nature 442, 1008 (2006). 
[12] E. Liang et al, Astrophys. J. 662, 1111 (2007); K. Murase et al, Astrophys. J. 651, L5 
(2006); N. Gupta and B. Zhang, AstroParticle. Phys. 27, 386 (2007) 



[13] S. Razzaque, P. Meszaros and E. Waxman, Phys. Rev. D 68, 083001 (2003) |astro-ph/0303505] . 
[14] S. Razzaque, P. Meszaros and E. Waxman, Mod. Phys. Lett. A 20, 2351 (2005) 

|astro-ph/0509729| . 
[15] S. Sahu and B. Zhang, Res. Astron. Astrophys. 10, 943 (2010) |arXiv:1007.4"582] [hep-ph]]. 
[16] V. D. Barger, K. Whisnant, S. Pakvasa and R. J. N. Phillips, Phys. Rev. D 22, 2718 (1980). 
[17] C. W. Kim and W. K. Sze, Phys. Rev. D 35, 1404 (1987). 
[18] H. W. Zaglauer and K. H. Schwarzer, Z. Phys. C 40, 273 (1988). 
[19] S. T. Petcov, Phys. Lett. B 191, 299 (1987). 
[20] H. Lehmann, P. Osland and T. T. Wu, Commun. Math. Phys. 219, 77 (2001) 

fhep-ph/0006213). 



[21] P. Osland and T. T. Wu, Phys. Rev. D 62, 013008 (2000) |h ep-ph/9912540] . 

[22] T. Ohlsson and H. Snellman, J. Math. Phys. 41, 2768 (2000) [Erratum-ibid. 42, 2345 (2001)] 



|hep-ph/ 9910546). 
[23] T. Ohlsson and H. Snellman, Phys. Lett. B 474, 153 (2000) |hep-ph/9912295] . 



[24] T. Ohlsson and H. Snellman, Eur. Phys. J. C 20, 507 (2001) |hep-ph/0103252] . 

[25] A. I. MacFadyen, S. E. Woosley and A. Heger, Astrophys. J. 550, 410 (2001) 



|astro-ph/9910034l. 



[26] C. D. Matzner and C. F. McKee, Astrophys. J. 510, 379 (1999) |astro-ph/99i0034] . 

[27] T. Shigeyama and K. Nomoto, C. D. Matzner and C. F. McKee, Astrophys. J. 360, 242 

(1990). 
[28] D. Arnett, Astrophys. J. 383, 295 (1991). 
[29] R. A. Chevalier and N. Soker, Astrophys. J. 341, 867 (1989). 
[30] O. Mena, L Mocioiu and S. Razzaque, Phys. Rev. D 75, 063003 (2007) |astro-ph/0612325j. 



19 



